Temperature of a Compressed Bubble with Application to Sonoluminescence 
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The rise in temperature from the adiabatic compression of a bubble is computed in thermodynamic 
mean field (van der Waals) theory. It is shown that the temperature rise is higher for the noble 
gas atoms than for more complex gas molecules. The adiabatic temperature rise is shown to be 
sufficient for producing sonoluminescence via the excited electronic states of the atoms. 
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I. INTRODUCTION 

There has been considerable interest in the dynamics 
of single bubble sonoluminescence ||, ^, ^. It is found 
that noble gases (within the bubble) greatly enhance the 
emitted light intensity when compared to bubbles filled 
with other gases |, |, 0, |, |, |l0|, |ll|. The theoretical 
interpretation of this experimental fact will be discussed 
in the work which follows. 

The fluid mechanics of bubble motion [|l2[ ^ 
is thought to play a central role in the dynamics of 
sonoluminescence. The starting point of theoretical 
bubble dynamics has been the Rayleigh-Plesset equa- 
tion [|l5[ |l6[ 0. An important feature of the exper- 
imental bubble fluid dynamics is the exceedingly high 
temperatures |l^, that can be achieved by merely 
decreasing the radius of a bubble. This is especially true 
when noble gas atoms are contained within the bubble 
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Although the dynamics of bubble 
compression appears to be intimately connected to the 
high noise temperature of emitted light pulses, the pro- 
cess of light emission has not yet been firmly established 

Our purpose is to calculate the maximum temperature 
of the contents inside the bubble on the basis of fiuid 
mechanics ^ |4| using a mean field (van der Waals) 
thermodynamic equation of state for adiabatic compres- 
sion. We find that the resulting temperatures are suf- 
ficiently high to produce light from sonoluminescence if 
not from mere incandescence. We also estimate the con- 
centration of atoms with excited electronic states which 
arise from an adiabatic compression of the fluid within 
the bubble. 

In Sec. II, we derive the thermodynamic identities re- 
quired to compute the temperature rise during adiabatic 
compression of a bubble. The mean field (van der Waals) 
thermodynamic equation of state is discussed in Sec. Ill 
wherein numerical results are presented for the case of ar- 
gon. Arguments are presented which explain why noble 
gas atoms are particularly suited for achieving high tem- 
peratures for small bubble radii. In Sec. IV, a Lagrangian 
formulation of the bubble dynamics will be employed to 
discuss in more detail the role of adiabatic processes. The 
internal chemistry|Q of the bubble will be discussed in 
Sec.V, while the physical aspects of the chemical reaction 



will be further explored in the concluding Sec. VI. 



II. THERMODYNAMICS 

If e(s, v) is the energy per particle of a fiuid inside the 
bubble as a function of volume per particle v and entropy 
per particle s, then 



de = Tds — pdv, 



(1) 



where T and p denote (respectively) the temperature and 
pressure of the fluid. Alternatively, one may employ the 
free energy per particle 



which obeys 



.f{T,v)^ini{e-Ts{e,v)}, 



df = —sdT — pdv. 



(2) 



(3) 



Condensed matter often (but not always) exhibits a 
lowering of the temperature when adiabatically expanded 
and a temperature rise when adiabatically compressed. 
The thermal coefficient describing such a process may 
defined as 



dT 

dv 



(4) 



From Eq.(|l|) it follows that 



dT 
dv 



dp 
ds 



dp 
df 



dT 
ds 



Eqs.(P and (^) imply 
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where 



Cv = T 



dp 
df 



ds_ 
df 



is the specific heat at constant volume. Similarly 



= T{^ 
' dT 



(6) 



(7) 



(8) 
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is the specific heat at constant pressure which obeys 

On the one hand, we have from Eqs.(|), (|) and (0) that 

(10) 

On the other hand, with 



dv }^ cA dT 



(11) 



Eqs.(^, die]) and (y) imply 

dT 
dv 
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(12) 



The temperature rise corresponding to an adiabatic com- 
pression is then calculated from 



(13) 



where (dT/dv)^ may be calculated from cither Eq.(p^) 
or Eq.(y). 

For the case of an ideal gas with Pqv = fc^T, 7 — const 
and TPo = 1, Eqs.(|l2l) and (P) imply 



T, 



T,, 



^ / {ideal gas) 



(7-1) 
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Compression of an Ideal Gas 




Y=(5/3) 
Y=(7/5) 



where A'' particles in a bubble of radius R are described 
by the volume 



V 3 



(15) 



The adiabatic compression of an ideal gas is plotted in 
FIG.|l|. For noble gas atoms (7 = 5/3) the rise in temper- 
ature is substantial If the radius drops by a factor 
of five, the ideal gas can exhibit a temperature rise from 
(say) room temperature to a temperature T ~ 7 x 10'^ K 



, ^ which is approximately the surface 
temperature of the sun. For an ideal gas of tumbling di- 
atomic molecules (7 = 7/5) the rise in temperature upon 
adiabatic compression is much less steep. 

The general situation follows from the thermodynamic 
Eqs.(p^ and (p^. For a complicated molecule, the inter- 
nal degrees of freedom of the constituent molecules give 
rise to a high specific heat c^. The adiabatic compres- 
sion of a fluid with a high heat capacity yields an only 
moderate rise in temperature. On the other hand, the 
noble gas atoms have at most internal electronic degrees 
of freedom and these lie dormant for reasonably low tem- 
peratures. Thus, noble gas atom fluids have a low heat 
capacity and a resulting high temperature rise on an adia- 
batic compression. These qualitative considerations give 
more than just a little insight into why noble gas atoms 
(in particular Ar atoms) have played such an important 
role in laboratory experiments on sonoluminescence. 

The sonoluminescent bubble can be so compressed that 
the fluid contained within the bubble is far from ideal. 
However, the above thermodynamic considerations are 
general and can be applied to fluids which are far from 
being ideal gases. In the mean field (van der Waals) the- 
ory of argon fluids, the equations of state can be worked 
out in detail. 



III. VAN DER WAALS THEORY 

The free energy per particle of an ideal gas of noble 
gas atoms may be written as 



fo{T,v)^-kBTlni^-^j, (16) 

where the thermal wavelength of an atom of mass m is 
given by 



At — 



mkBT 



(17) 



FIG. 1: For the adiabatic compression of an ideal gas in a 
bubble from radius Ri Rf < Ri, the temperature rise is 
described by Ti Tj > Ti. For the case of noble gas atoms 
(7 = 5/3) the temperature rise is much more steep than the 
case of tumbling diatomic molecules (7 = 7/5). 



In mean fleld theory, the hard core of the two atom 
potential is included by the "free volume" replacement 
V ^ V — h while the attractive part of the potential is 
simulated by subtracting (from the free energy) an at- 
tractive mean potential term {a/v) proportional to the 
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fluid density. Thus, the "mean field" theory modifica- 
tion to the ideal gas Eq.(^6|) is given by 



'ksTln 



V ^ b 



(18) 



The mean field theory pressure p = —{df/dv)T which 
follows from Eq. (n8) is that of van der Waals 



P 



(19) 



while the specific heat at constant volume follows from 
Eqs.(0) and (|l|) as 



_9£ 

Eqs.(^, (|l|) and (|^ imply that 
'dT\ ( 2T 



(20) 



dv 



3{v - b) 



(van der Waals). (21) 



Finally, Eqs.(|l3D and (21) yield the mean field theory 
adiabatic temperature rise in a noble gas atom fluid upon 
adiabatic compression; It is 



(van der Waals) 



2/3 



Vf - b, 
(i?//i?,)3-y3 



Adiabatic Temperature 



2/3 



,(22) 
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FIG. 2; In the adiabatic compression of a fluid within a spher- 
ical bubble from radius Ri to radius Rf, the temperature rise 
from Ti to T/ is plotted. The temperature rise from room 
temperature when interactions are included (as in the van 
der Waals theory) is large compared with the ideal gas case 
but only if the compression is substantial. 



where Eq.(|lJ) has been invoked and the initial volume 
per particle Vi determines 



y = (b/v,)'/^ = i^/A) 



(23) 



In Eq.(|2|), 6 = $3 ^ ^3^ Pqj. j^j.goj^ [|7|, C « 

3.7664 A and for an initial state at room temperature 
and pressure A « 34.356 1. Thus y{Ar) « 0.1096. A 
plot of Eq.(|2|) is shown in FIG.|. 

If the final bubble radius i?/ is only (say) « 10% of the 
initial radius Ri Q, ^ then in the van der Waals 
model the temperature rise can be enormous (say) from 
room temperature to a final temperature of « 7.5 x 10* K 
at which point the fluid would appear to be incandescent. 



IV. LAGRANGIAN DYNAMICS 

For the purpose of better understanding the role of 
adiabatic processes in bubble dynamics, we here con- 
sider such dynamics from an analytical mechanics view- 
point. The fluid mechanics of a spherical bubble can be 
described by a Lagrangian of the form 
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L{R,R) = ^M{R)R^ -U{R) ~Uext{R) 



(24) 



where the external driving potential for the bubble dy- 
namics is denoted by lAext and a Rayleigh dissipation 
function Wn T is employed having the form 



2T{R,R) = T{R)R^. 
The dynamical equation of motion reads 



d_ /dL 

dt VdR 



dL 
OR 



dR 



(25) 



(26) 



The various terms contributing to Eqs.(p4|) and ( P5[) will 
now be explored. 

During an expansion or compression of a spherical bub- 
ble, the fluid velocity field outside the bubble is 



R^R for (|r| > R) (27) 
while the fluid velocity field inside the bubble is 

r for (|r| < R). (28) 

The (presumed) incompressible flow outside the bubble 
requires a kinetic energy 



\r\>R 



R'R 



2TTpR^if, 



(29) 
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while the flow of the fluid inside the bubble requires a 
kinetic energy 



M 



2(47ri?3/3) 



W?d?v 



3M 
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rR 
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Jo 


R 



\<R 
2 



(47rr^dr) 



= {3/W)MR^, 



(30) 



where M denotes the mass of fluid (presumed) conserved 
within the bubble. From Eqs.(|2^) and ( pO| ) one finds a 
total fluid flow kinetic energy of 



for which 



M{R) = -M + At:pR^. 
5 



(32) 



Eq. (p2[) establishes the effective mass M{R) to be used 
in the Lagrangian Eq.(p^. 

For the potential terms in the Lagrangian Eq. ( p4[ ) , one 
may write a sum of the internal energy of the fluid within 
the bubble and a term proportional to the surface area 



47rcri?^ 



(33) 



The external potential may be written in terms of the 
fluid pressure P exerted on the bubble by the external 
fluid; i.e. 



where 



r(i?) - 1671 ( 77 + ) i?. 



(38) 



Eq.(p8[) establishes the friction coefficient T{R) to be used 
in the dissipation function Eq. (|25|). 

Employing Eqs.(0), (|4[]2|), (p[||) and ^ we find 
the following equation of motion: 



{l + ai)RR 




(31) where 



(39) 



(40) 



is the kinematic viscosity of the liquid outside the bubble, 
and ai for i = 1,2 in Eq.(39) represent the corrections 
to the usual Rayleigh-Plesset equation ||l^, ^ |l^ due to 
motions within the bubble. In detail, 



(41) 



In the usual approximation, a2 ~ 0, all of the dissipative 
heating is from outside the bubble, and the compressions 
and/or expansions of the bubble are strictly adiabatic. 
For the internal heating of fluid inside the bubble, a2 > 0, 
the bulk viscosity ( of the fluid inside the bubble will give 
rise to additional internal heating. 



l^ext {R) 



UttR^P 



(34) 



Eqs.(|2l), (H) and (||) completely specify the La- 
grangian in Eq. (|2^) . 

To compute the dissipation function, wc note for the 
incompressible flow (with viscosity 77) [^Sj outside the 
bubble 



\r\>R 



l6nriRR^, 



(35) 



where Eq.(g7|) has been invoked. For the flow inside the 
bubble, only the bulk viscosity C |Q enters into the com- 
putation, i.e. 



2T 



C / (divvfcfr 

J\r\<R 

l2nCRR^, 



where Eq.(p^) has been invoked. From Eqs. 
one finds a total dissipation function 

2^ = 2^„ + 2^o„t = r(i?)i?2 



(36) 
|) and ^) 

(37) 



V. ELECTRONIC EXCITATION OF ATOMS 

Let us now further consider the case in which the fluid 
inside the bubble is made up of argon atoms. We wish 
to consider the physical chemistry of a simple reaction in 
which the electron states of argon atoms are excited; i.e. 



Ar* 



(42) 



where Ar represents an atom in the ground electronic 
state while Ar* represents an atom the atom in an excited 
electronic state. The concentration of atoms in an excited 
electronic state is defined as 



Nav 



Nat + Nav 



(43) 



The second law of thermodynamics dictates that chem- 
ical equilibrium arises when the energy reaches a mini- 
mum for constant entropy; i.e. 



inf e{s,v,x). 

0<x<l 



(44) 



When the fluid in the bubble is out of chemical equilib- 
rium, Eq.(|^) is generalized to 



de — Tds — pdv + jj-dx, 



(45) 
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where 



VI. CONCLUSIONS 



(46) 



If a; << 1, then it is sufhciently accurate to use the con- 
ventional dilute solution theory wherein 



fl — (p + ksTlnx 



(47) 



and the activation free energy (j) is independent of x. 
From Eqs.(44) and (Ea), it follows for chemical equilib- 
rium that fL — 0. Eq.(p1) then yields the thermal activa- 
tion formula for concentration 



(chemical equilibrium) 



(48) 



In thermodynamic mean field theory, the activation 
free energy may be written in terms of the energy states 
of a free argon atom as 



= —knThi 



E 



E>Eo 



Wee 



-E/kBT' 



(49) 



where we the degeneracy of the electronic energy level 
E and Eq is the electronic ground state energy. For argon 
it is known|^ that the first two electronic excited states 
are given by 



E 



/2 



)4s, J = 1 -Eo 



11.63 eV 
11.82 eV. (50) 



Including the above electronic states in the computation 
of the activation free energy in Eq.(^) and employing 
Eqs.(^8|) and (^) we plot the concentration of excited 
states as a function of bubble radius in FIG.|^. 
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FIG. 3: As bubble containing pure argon is adiabatically com- 
pressed from an initial radius Ri at room temperature and 
pressure to a final radius Rf < Ri, the concentration x of 
electronically excited argon atoms sharply increases. 



It has been shown that the adiabatic compression of 
noble gas atoms gives rise to surprisingly high tempera- 
tures. The statistical thermodynamics of the processes 
may viewed as follows: (i) The entropy of the fluid inside 
the bubble is proportional to the logarithm of the num- 
ber of states, (ii) When the bubble is compressed, due 
to the hard core part of the inter-atomic potential, there 
is an excluded spatial volume. The compression then 
eliminates much available spatial volume part of phase 
space which would in principle lower the entropy, (iii) 
To keep the entropy constant during the adiabatic com- 
pression, the momentum part of phase space must ex- 
pand as the position part of phase space contracts. Thus 
there will exist an increased kinetic energy. A high mean 
kinetic energy (in classical statistical thermodynamics) 
is equivalent to a high temperature. This explains why 
an adiabatic bubble compression can lead to such a high 
temperature. 

That noble gas atoms achieve even higher adiabatic 
compression temperatures than the temperatures which 
would be achieved with more complex molecular sub- 
stances has also been explained. The compression lowers 
the spatial parts of phase space so other parts of phase 
space must be opened up to keep the entropy constant 
in an adiabatic process. For a fluid of noble gas atoms, 
momentum phase space is opened up giving rise to high 
kinetic energy and high temperature. If the fluid inside 
the bubble contains molecules with structure, then other 
parts of phase space can open up to make the lost spa- 
tial part of phase space. These include internal tumbling 
motions of molecules as well as internal vibrational ex- 
citations. Not all of the required new regions of phase 
space involve a rise in center of mass kinetic energy of 
the molecule which is equivalent to a rise in tempera- 
ture. Thus, a fluid composed of complex molecules will 
not exhibit as high a compression induced temperature 
rise as will a fluid made up of noble gas atoms. From a 
mathematical point of view, the specific heat Cy in Eq. (|^) 
describes how much phase space is opened up by a given 
temperature rise. Each of the internal degrees of freedom 
in a complex molecule raises c^. The higher the specific 
heat, the lower the temperature rise in Eqs.(|0l) and (|l|). 

The bubble dynamics has been described through a 
Lagrangian, including the Rayleigh dissipation function 
yielding the usual results apart from small correction 
terms arising from fluid flow within the bubble. Since 
the energy for fixed entropy enter into the Lagrangian as 
potential energy terms, the adiabatic compression holds 
to a high degree of accuracy. 

The concentration of electronically excited atoms at 
compression induced high temperatures can be approxi- 
mately estimated in mean field theory. One uses an elec- 
tronic partition function which is strictly valid only for 
ideal gases. When the fluid is compressed, there will in 
reality be a shift in electronic energy. The excited elec- 
tron shells will be somewhat delocalized to neighboring 
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atoms. In fact, one might well have coherent electron garding the coherence of sonoluminescent optical pulses 
states in such a system giving rise to superradiance dur- is still unclear and requires further elucidation, 
ing the light emission pulse. However, the situation re- 
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